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Recently, a method has been proposed to obtain accurate predictions for low-temperature prop- 
erties of lattice spin glasses that is practical even above the upper critical dimension, dc — 6. This 
method is based on the observation that bond-dilution enables the numerical treatment of larger 
lattices, and that the subsequent combination of such data at various bond densities into a finite-size 
scaling Ansatz produces more robust scaling behavior. In the present study we test the potential of 
such a procedure, in particular, to obtain the stiffness exponent for the hierarchical Migdal-Kadanoff 
lattice. Critical exponents for this model are known with great accuracy and any simulations can 
be executed to very large lattice sizes at almost any bond density, effecting a insightful comparison 
that highlights the advantages - as well as the weaknesses - of this method. These insights are 
applied to the Edwards- Anderson model in d = 3 with Gaussian bonds. 
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I. INTRODUCTION 

The exploration of low-temperature properties of dis- 
ordered systems remains an important and challeng- 
ing problem Systems in this class possess a low- 
temperature glassy state with a glass transition at some 
temperature Tg > 0. They are characterized by a com- 
plex (free-)energy landscape in configuration space 0,S 
U with a hierarchy of valleys and barriers whose multi- 
modal structure impedes the progression of any dynam- 
ics towards equilibration, causing tantalizing phenom- 
ena, such as trapping and jamming on intermediate time 
scales, and aging on long time scales. An understand- 
ing of such systems is of paramount importance as these 
phenomena are observed for a large class of materials as 
well as for biological systems 0- 

The paradigmatic model for the study of such phenom- 
ena is the Ising spin glass, either on a finite-dimensional 
lattice (Edwards- Anderson model, EA '3]) or some other 
random network structure. Disorder effects arise typi- 
cally via quenched random bonds, random local fields, 
or merely the randomness of the network itself, each can 
lead to conflicting constraints which leave variables frus- 
trated. It is believed that a proper understanding of 
static and dynamic features of EA may aid a description 
of the unifying principles expressed in the wider class of 
realistic problems 0- 

Unfortunately, after 30 years of research there is still 
not consensus on whether the subtle mean-field picture 
for Ising spin glasses derived long ago H, 0, 0| has 
any relevance for real-world materials at low temper- 
atures 0. Progress towards such an understanding 
is slow due to the peculiar structure of the problem 
at hand. The intricate multi-modal low-energy land- 
scape puts the computational effort needed to deter- 
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mine thermodynamic observables usually into the class of 
NP-hard combinatorial optimization problems 12\ know 
from computer science, for which worst case computa- 
tional costs increase faster than any power of the system 
size. (Interestingly, it is in this area of computer science 
itself where the mean-field theory at T = has had a 
most significant impact so far ji^]-) 

Most insights into finite-dimensional systems has 
thus been gained through alternative computational ap- 
proaches to elucidate low-temperature properties. Aside 
from methods designed to expedite some thermodynami- 
cally correct algorithm, such as parallel tempering 
or the waiting time method [l5j, one focus area has 
been the use of optimization heuristics directed toward 
fully enumerati ng grou nd-state configurations by any 
means |1 H [H [lllll . These ground states should 
provide the basis for the thermodynamic behavior of the 
system at and near T = 0. But due to the NP-hardness, 
even heuristic methods become unreliable when system 
sizes exceed about « 10^ variables, often too small to 
draw safe conclusions or to sufficiently discriminate be- 
tween theoretical ideas ll^ • 

This fact is illustrated by the determination of the stiff- 
ness exponent, often labeled y or 6 (]i^23j, a fundamental 
quantity assessing low-temperature energy ffuctuations: 
a positive value of y, as found in EA for d > 3, de- 
notes the increase in the energetic cost (i. e. "stiffness") 
accompanying a growing number of variables perturbed 
from their position in the ground state. The rise in en- 
ergetic penalty paid for stronger disturbances signals the 
presence of an ordered state. In turn, for systems with 
y < such order is destabilized by arbitrarily small fluc- 
tuations. 

In this paper, we will extract the stiffness exponent 
from the response induced through defect-interfaces ^3 ■ 
These can be created by fixing the spins along the two 
faces of an open boundary in one lattice direction. The 
ground state configuration with energy Eq of an in- 
stance is first determined for a random fixing of those 
boundary spins, then the energy Eq is obtained for the 



same instance and the same fixing, but with all spins 
reversed on one of the faces. Hence, the interface en- 
ergy Ai? = Eq — Eq created by the perturbation on the 
boundary is sampled, and its distribution P(AE) deter- 
mined. If a system is glassy, the typical energy scale 
involved, here re presented by the width of the distribu- 
tion, a{AE) = y/{AE'^) - (AE)'^, should grow with the 
size of the perturbation, say, the linear extend of the 
boundary, L, as 



a{AE) - i» 



(1) 



Accurate determination of this exponent in d = 3 has 
long been elusi ve, with values given between y w 0.19 0, 
1131 to « 0.27 m. While it was save to say that y > 0, 
its value could be at best given as j/3 = 0.2(1). In terms 
of sorting out theoretical models, this is too inaccurate 
to draw solid conclusions [2^. 

This limited accuracy originates with trying to fit data 
over too small a range of system sizes L. Scaling in 
this range is further beset by slowly decaying correc- 
tions in L which make it hard to decide even where scal- 
ing sets in [2^ [S^. The remedy proposed then is to 
possibly extend scaling by considering bond- diluted lat- 
tices [23, [2^. As has been argued already in Ref. [s^, 
as soon as the percolation window for bond densities 
p > Pc has been exceeded at some lattice size, the long- 
range properties of the giant component - and hence, of 
the spin glass defined on it - are essentially compact, 
representing the T — fixed-point of the fully connected 
lattice. This implies that there the exponent y is inde- 
pendent of the bond density p. These features of bond- 
diluted glasses are summarized in the phase diagram in 
Fig. n The independence of y from p appears to hold 
also for Gaussian bonds in d = 2 where y < (no 
spin glass phase) and universality with respect to bond 
disorder is violated pllsills^. 

While we focus on the defect energy at T = here, 
bond-dilution may also be an effective means to study 
other observables | i34l] . For determining T — proper- 
ties, this approach makes the treatment of larger lattices 
sizes practical in as much as exact algorithms can be 
devised that allow the elimination of a large number of 
variables, whose state becomes entrained to other vari- 
ables in a predictable way in the ground state . In this 
way, systems with more than 10^ variables have been "re- 
duced" to remainder graphs consisting of no more than 
a few hundred variables that are amenable to standard 
optimization methods |23] . Clearly, though, due to those 
almost trivial variables, the information contend about 
the asymptotic behavior captured by such graphs may 
be limited, and a price has to be paid through extended 
transient behavior. As we hope to demonstrate here, in 
the end much is gained in this method, although at the 
extremely large lattice sizes that can be reached with 
the Migdal-KadanofF approximation used here, dimin- 
ished returns are obtained because of a lack of knowledge 
on how scaling corrections vary with bond density. 

As we will recount in Sec. lIIII Refs. |3^,32i| suggest to 
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FIG. 1: Phase diagram for bond-diluted spin glasses at low 
temperatures. Lattice spin glasses of dimension d > 3 possess 
a spin glass phase (SG) below the glass transition temperature 
Tg > 0. This phase should persist even for diluted lattices, as 
long as the percolating cluster of connected spins is compact, 
i.e. for bond densities p sufficiently above Pc |3^ . 
Below Pc, the collection of small clusters can only respond 
paramagnetically (PM) on long lengthscales for any T > 0. 
Depending on the details of the bond distribution, the T = 
transition occurs at p* = Pc (continuous bonds) or at p* > pc 
(discrete bonds), respectively determining the slope of the 
phase boundary Tg{p) at p* . 



generalize Eq. ^ to 



L \ 



for L » 1 and ^(p) ~ (p - p*) > 1. The scahng 
fmiction / was chosen to be constant for L £,{p)- The 
limit p ^ p* towards the T = 0-transition between the 
spin glass and paramagnetic regime (see Fig. ^ is inter- 
esting in its own right, and will be investigated further 
in Ref. js^l- Yet, to elucidate properties of the glassy 
regime, the "window of opportunity" for our method ap- 
pears to be at intermediate bond densities: p has to 
be sufficiently smaller than unity for our reduction al- 
gorithm to be efficient, but also sufficiently above p* to 
attain system sizes L 3> S,ip)- Analyzing our numeri- 
cal data here suggests that in this window the condition 
S.ip) ~ ip-p*)~'^' > 1 assumed in Eq. does not hold, 
making ^ (p) a more general function of p with unknown 
corrections to its singular part. 

Operationally, we thus propose a naive Ansatz for a 
collapse of the numerical data valid for the asymptotic 
regime L ^ £^{p) 1 only. 



(7{AE) = f{oo) xy 



x = L{p-p*r 



(3) 



Since the numerically accessible data appears to violate 
^(p) 3> 1, the parameters p* and v* here merely facilitate 
the data collapse by fitting a more general function ^(p), 
and can not be expected to yield accurate predictions for 
the critical values in Eq. Q . The successful implementa- 
tion of this Ansatz for the Migdal-Kadanoff hierarchical 



lattice here lends credibility to the findings for yd of the 
EA in Refs. (HHI. 

Finally, it is remarkable that this approach works sub- 
stantially better for a discrete ± J bond distribution than 
for Gaussian bonds, for which scaling corrections due to 
a small ^(p), and even due to finitc-L transients, are far 
more significant. Just as for the ±J data, the similar- 
ity between our Migdal-Kadanoff and our EA data with 
Gaussian bonds presented here is striking. While this is- 
sue eventually deserves more thorough investigation, we 
can speculate on the origin of these strong transients. 
One may remember that in d = 1, a{AE) ~ L^^ for 
continuously distributed bonds with finite P(0), while 
a{AE) ~ L'^ for ±J bonds. The interface may set- 
tle on the extremely weakest within a more widely dis- 
tributed set of continuous bonds, while for ±J it must 
break a bond of order unity. Similarly, on short ranges in 
higher dimensional lattices, particularly dilute ones, with 
a Gaussian distribution peaked at J = 0, the violation of 
heavy bonds can be deferred initially that an extended 
defect eventually demands on larger scales L. This naive 
argument is supported by the fact that at equal p and 
L, defect energies a{AE) are typically somewhat smaller 
for Gaussian than for ± J bonds, as will be seen in Fig. 13 

In the next section we introduce the Migdal-Kadanoff 
hierarchical lattice, followed in Sec. IIIII by a discussion 
of the scaling arguments leading to Eq. Q and to our 
Ansatz in Eq. (|3J) used to obtain the stiffness exponent 
on bond-diluted lattices. Sec. IIVI contains our discussion 
of the Migdal-Kadanoff data, followed by the application 
of the gained insights to the diluted EA with Gaussian 
bonds in Sec. Finally, Sec. IVII contains our conclu- 
sions. 



II. MIGDAL-KADANOFF HIERARCHICAL 
LATTICE 

To test the scaling Ansatz in Eq. Q for the stiff- 
ness exponent, we consider here the bond-diluted hi- 
erarchical lattice (see Fig. EJ, obtained in the Migdal- 
Kadanoff bond-moving scheme 38] for low-dimensional 
spin glasses. These lattices have a simple recursive yet 
geometric structure and are well-studied [23, ISJ, [S^, Eol ■ 
Most importantly, their ground states can be obtained in 
polynomial time at any bond density, and we can discuss 
our results independent of any systematic bias introduce 
by a subsequent optimization process that may be re- 
quired for more complicated models like EA 'ss'l . A most 
interesting property of these lattices is the curious fact 
that the scaling of its defect energy distribution, giving 
rise to the stiffness exponent, behaves in all respects very 
similar to that measured for actual three-dimensional lat- 
tices USES- 

To generate a hierarchical lattice, starting from gen- 
eration 7 = with a single link, at each subsequent 
generation / all links from 7—1 are replace with a 
new subgraph, as described in Fig. [21 The structure 
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FIG. 2: Diagram for the recursive algorithm to calculate spin 
glasses on hierarchical lattices. Drawing bonds J/_i randomly 
from a sufficiently large pool P/_i(J) at generation J — 1 
(top), a lattice of generation I is formed, then bond-diluted 
to density p, and finally "reduced" 35] to an effective bond 
Ji between the root spins (squares) , which is add to the pool 
Pi (J) (bottom). 

of the subgraph arises from the Migdal-Kadanoff bond- 
moving scheme in d dimensions, and has 2"^ = 8 links 
for d = 3 here. Thus, a hierarchical lattice of gener- 
ation / has {2'^)^ links (if undiluted), thus correspond- 
ing to a d-dimensional lattice of "length" L = 2^ and 
n = 2-h2'^-i(L''-l)/(2'*-l) = OiL"^) vertices. While the 
average connectivity is 4 — 2'^^'^, the two root-vertices 
from generation 7 = themselves obtain in generation 7 
a connectivity of ^ 2*^'*~^^'^-^~^^ In turn, 2"^^^^ vertices, 
7 in 8 for d = 3, are only two-connected. 

The diluted hierarchical lattice percolates when there 
is a path between the two root- vertices, representing the 
boundaries of the system. This notion leads to a simple 
recursion relation for the percolation threshold by count- 
ing the weights of all diluted subgraphs from Fig. [3 that 
percolate, i. e. connect right and left vertex. In d = 3 
one gets 

Pi+i=Apj-6p^j + 4p^j-pl (4) 

which has a non-trivial stationary point at pc — 
0.2818376366. There can not be long-range correlated 
behavior, such as spin glass orde ring , for any p < pc- 
It has been pointed out by Ref. [33 that a spin glass 
on a hierarchical lattice with ±J-bonds exhibits a crit- 
ical transition between a paramagnetic and a spin glass 
phase for bond densities at p* = 0.31032. While below 
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Pc disconnected clusters clearly prevail and prevent long- 
range correlations, even for < p < p* such correlations 
remain suppressed due to the cooperative behavior in the 
bond structure pervasive in the lattice, leading to cancel- 
lations that additionally disconnects subgraphs at some 
higher level of the hierarchy. In contrast, for a continuous 
bond distribution, such as the Gaussian bonds discussed 
below, any such cancellations would be unlikely, leading 
in this case to p* — Pc- 

III. DEFECT-ENERGY SCALING FOR p ^ p* 

Following the discussion in Refs. [ssL Is^ . for bond- 
diluted lattices at p — > p* we have to generalize the scal- 
ing relation in Eq. |^ for the defect energy to 

a{AE)L,p ^ y{p)Lyf (^^) , (5) 

where y ^ {p " P*Y is an effective surface tension and 
^(p) ~ (p — p*) '^ is the correlation length for the cross- 
over into glassy behavior. Note that Eq. lO requires 
both, L and ^, to be large compared with the unit lattice 
spacing to avoid further scaling corrections. The scaling 
function f{x) is defined to be constant for large argu- 
ment, L > C > 1. 

For ^ ^ L '3> 1, Eq. ^ requires that f{x) ~ x^^ for 
X — * to satisfy either the vanishing of a with L in case 
of p* = Pc (Gaussian bonds), or its scale invariance at a 
p* > Pc (±J bonds). Clearly, due to the tenuous fractal 
nature of the percolating cluster at p* = Pc, no long- 
range order can be sustained, defects possess a vanishing 
interface, and one may expect that 

a{AE)L,r --Ly-, (6) 

where yp <0 js^ . To cancel the p-dependence at p = p* , 
Eq. ||SJ) requires y+fi = yp and t+fiu = 0, i. e. t = iyy+4>, 
setting (/) = —vyp [s^. As a result, we obtain Eq. 
In contrast, for p* > Pc, the percolating cluster appears 
compact on scales L ^ (p* —pc)~^, and a can not vanish 
for further increasing L. Yet, neither can a increase at 
p* , by definition. Hence, a remains scale invariant, and 
correspondingly cj) and yp vanish in Eqs. (0) and © for 
± J bonds. 

In the case of continuous bond-distributions P{J) with 
finite P(0), where p* = Pc, spin glass ordering can only be 
noticed on scales at least as large as the correlation length 
^ ^ {p — Pc)~'' associated with percolation, suggesting 
u* = V. In case of ±J-bonds, the cooperative cancel- 
lations between discrete bonds mentioned above further 
weaken order, leading to v* > v for the T = glass tran- 
sition at p* . Accordingly, Ref. 33] finds v = X.TXIA, and 
V* = 1.5373 for the hierarchical lattice in d — i. 

Finally, at the cross-over ^ ^ L, where the range L 
of the energy excitations (t{AE) reaches the percolation 
length and spin glass order ensues, Eq. yields 

a{AE)^^,y, ^{p- p*f e/(l) -{P- P*f • (7) 



One can associate a characteristic temperature with this 
cross-over by (3a{AE)^(^p-^^p ~ 1; for temperatures above 
this T = 1//3, thermal fluctuation destroy spin glass or- 
der. This suggests a relation between bond density and 
the glass transition temperature: 

T,(p)^(p-p*)^ (8) 

Eq. © defines as the "thermal-percolative cross-over 
exponent" ,36j, which specifies the details of the phase 
boundary near pc (or p*) in Fig. ^ For the continuous 
distribution, where p* — Pc, 4> would be in general non- 
trivial, while for ±J at p* > pc it appears that (p = 0, 
indicating a jump in the phase boundary Tg{p) at p* . 

The exponent (p could be of importance, since it may 
be experimentally accessible 0, li^ while the relation 
(j) = —ypv provides a simple computational determina- 
tion in terms of Eq. © and the well-known percolation 
exponents v. This connection will be considered in a 
forthcoming publication [3^ . 

In the following, we will explore some of these rela- 
tions numerically for the hierarchical lattice, for which 
large sizes L can be obtained. But the central purpose 
of this paper is to probe Eq. ©, which has been used 
in Refs. [23, ll^ to provide accurate predictions for the 
stiffness exponents y^, fundamental for describing low- 
temperature excitations in spin glasses. 

IV. NUMERICAL RESULTS FOR 
HIERARCHICAL LATTICES 

Our numerical studies on the hierarchical lattice have 
been conducted with the algorithm described in Ref. . 
It is based on the evolution of bond-pools (see Fig. O of 
size Aj from generation / to generation / -|- 1, similar 
to the procedure already used in Refs. 0| . In our 
algorithm, though, an existing bond at generation / is 
replaced with a new bond (to keep Ai constant) for ev- 
ery k new bonds that are added at generation / —I. This 
procedure is legitimate in principle, as even neighboring 
sub-elements in the graph act independently and only ef- 
fect each other in a collective sense just as represented 
by the bonds replacing them in the next generation. The 
danger is that the diversity in the pool of those bonds 
is insufficient, leading to creeping spurious correlations 
and difficult-to-perceive drifts away from the true values 
of observables. Since the size of the pool of bonds has 
only a minor effect the on computational effort for fc = 1, 
we can re-run the same calculation repeatedly to a high 
number of generations using ever-larger pool sizes Ai un- 
til the data becomes insensitive to Ai. We have used this 
procedure to generate most of the plots in Figs. 13 

Further problems arise when we use this algorithm near 
the critical point p* . There, the bond distribution is torn 
between the trivial fix point of vanishing width in the 
paramagnetic phase and the true glassy state with subtle, 
long-rang correlations between bonds of the previous and 
following generations, and small fluctuations can severely 
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FIG. 3: Plot of the width cr( Ai?) of the defect energy distribu- 
tion as a function of systems size L for various bond fractions 
p using discrete bonds (top) and Gaussian bonds (bottom) in 
the Migdal-Kadanoff lattice. In each plot, dashed lines from 
bottom to top refer to p = 0.1, 0.2, . . . , 1. Solid lines refers 
to data with p = 0.28, 0.29, . . . ,0.34 on top, and p = 0.28, 
0.2818, and 0.29 on the bottom. In both plots, the data 
points with circles were obtained using the = 1 implemen- 
tation from Ref. [ssf . those with diamonds using fc = 2, and 
those with crosses near the respective p* required the exact 
algorithm. Note that at p* « 0.310 > pc on top, the slope 
is near vanishing, while for p* = Pc = 0.2818 on the bottom, 
the slope indicates a well-pronounced power-law decay. 



bias the evolution. At that point, we have to resort to 
slower regeneration rates with fc > 1, or even the exact 
algorithm for which fc = 2''. 

Our numerical studies, as shown in Fig.|21 confirm the 
picture described in Sec. IIIII To induce an interface in 
the MK as described in the Introduction we consider only 
the left-most and the right-most spin in each graph as 
entire boundary. This is meaningful, since each of these 
spins connects to 0(V) other spins, see Sec. Then, 
the defect energy of a lattice at generation / — 1 (size 
L — 2^"^) is simply (twice) the value of the bond of 
generation / replacing that graph (see Fig. ^ . In this 
sense, the defect energy a{AE) can be interpreted as 
an effective coupling between both sides of the defect 
interface |36j |: if that coupling strengthens with distance 
L, the system is in an ordered state, and vice versa. 



Note that in Fig. 13 for p < p* the data evolves towards 
the p = fix-point with (t{AE) — 0, while for p > p* 
it invariably evolves to the p = 1 fix-point with scal- 
ing behavior a{AE) ~ L^. Fig. O clearly suggests that 
p* = Pc = 0.2818 for the continuous Gaussian bond dis- 
tribution, while the discrete ± J bond distribution favors 
p* = 0.3103 > Pc- At p*, a varies distinctively different 
with L than even for bond densities p quite near to p* . 
In particular, the behavior of a at the respective p* ap- 
pears to confirm the predictions for the scaling in Eq. ® 
with Up = (i. e. = 0) for discrete bonds, and a non- 
trivial exponent we measure to be about yp — —0.9(1) 
for Gaussian bonds. 

In the following, we want to test the predictive power 
of the finite-size scaling Ansatz in Eq. Q as argued on 
the basis of Eq. (|2Jl. This Ansatz has been applied to 
the numerical measurements of the glassy state (p > p*) 
on large, bond-diluted lattices as recently proposed in 
Refs. ll^llSiJ. In this scheme, the properties of the zero- 
temperature fix-point are determined for a number of 
intermediate bond-densities p - sufficiently above p* to 
be glassy and sufficiently small to achieve large system 
sizes L. In particular, Refs. [23, focused on the de- 
fect energy for finite dimensional lattices. This procedure 
should be applicable also to other observables. 

First, we consider ±J bond-distributions, which had 
been used exclusively in Refs. 113, H^]. In Fig.0|we plot 
without any scaling all data we have obtained with the 
method above for p = 0.35, . . . , 1. It's worth mention- 
ing two features of that data: (1) scaling corrections for 
smaller L appear to change sign at around p « 0.45, 
suggesting that those corrections are weakest at interme- 
diate values of p instead of at p = 1, against expectation. 
(2) The data for p ^ 1 initially narrows (on this loga- 
rithmic scale) for equal increments (Ap — 0.1), but then 
exhibits an increased gap in the jump from p — 0.9 to 
p — 1. Both of these features have also been observed for 
the EA data in d = 3 [l^, HI]. The first feature seems 
insignificant here for data extending out to L > 10^, but 
the suppression of scaling corrections becomes extremely 
helpful when the maximal attainable L is small, as for the 
EA model. This fact has also been exploited in Ref. [33| . 
The second feature could be explained with the require- 
ment for the scaling Ansatz of ^ ~ (p — p*)""^ ^1, 
which may not be satisfied for any p too far from p* . In 
Ref. [3^, similar rapid variations in an observable (the 
overlap) for p — >• 1 have been observed, although a con- 
nection to the variations in the amplitude of a here is not 
clear. 

For this discrete bond-distributions, the scaling Ansatz 
in Eq. (0) simplifies, since we can assume yp = 0, see 
top of Fig. 131 Hence, we fit the data obtained for vari- 
ous bond densities p > p* hi Fig. ^ directly to the form 
proposed in Eq. where the scaling variable x is ad- 
justed to provide the best data collapse. The unknown 
scaling function f{x) has been replace by a fitting con- 
stant, /(oo), to capture its leading asymptotic behavior 
for large L. 
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FIG. 4: Plot of the raw data for a{AE) as a function of 
systems size L for bond fractions p* < p < 1 using discrete 
±J-bonds. Some of this data is already shown in Fig. |3 
Of note is that there appears to be no variation with p in 
the asymptotic scaling, that scaling corrections at small L 
are least noticeable for intermediate p, and that there is an 
anomalous gap between data for p = 0.9 and p = 1. 




0.01 0.1 1 10 100 , 1000 10000 100000 

L(p-py 

FIG. 5: Collapse of the data from Fig. |1] according to the 
scaling Ansatz in Eq. but using the exactly-known values 
for p* , u* , and y — 0.2552. Here, only /(oo), an overall 
amplitude, remains as a free fitting parameter. 



The advantage of the hierarchical lattice is that we 
already know the expected values for the parameters in- 
volved in such a fit i. e. p* = 0.31032, ly* = 1.5373 HI, 
and y w 0.2552 |33,|43|- Thus, first we can attempt to 
utilize this knowledge to collapse the data in Fig. 0] by 
fixing p*, V* , and y, only leaving /(oo) to be fitted. As 
can be seen in Fig. the data collapses exceedingly well 
even for small x, as long as p is not too large. Scaling en- 
sues quickly for a; > 1 , which would allow for an excellent 
fit, if one were to exclude data for p > 0.7. Remarkable is 
the quality of the collapse for data reaching below x < 1. 
It reveals a concave shape for f{x), which explains the 
finite-size corrections in Fig. 0] Data with lower p first 
rises slowly, than more rapidly for increasing L, before 
scaling settles in. In turn, data with higher p immedi- 
ately rises rapidly before settling into a slower asymptotic 
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FIG. 6: Collapse of the data from Fig. |1| according to the 
scaling Ansatz in Eq. Here, the values for f(x), p* , v* , 
and y are all determined from the fit of data in the scaling 
regime only, although all data is displayed. 



growth. Data with p near unity rises even more rapidly 
to overshoot the collapse in the scaling regime. 

Clearly, these observations about the intricacies of f{x) 
are far too subtle to be of any use in the more typical sit- 
uations where we have no prior knowledge of the param- 
eters. Worse yet, even with that knowledge, corrections 
due to the finite size of ^ = (p — p*)"" in particular 
appear to prevent a collapse of the data in much of the 
asymptotic scaling regime. Yet, it seems obvious that the 
data for all p > p* and large enough L exhibits scaling 
according to Eq. providing information we desire to 
exploit. The precise value of the parameters p* and v* 
are important for the behavior of f{x) near a; ~ 1, but 
have little effect on the scaling behavior for x 3> 1. Thus, 
instead of correcting the fit to extract accurate values for 
all parameters, we cut data that does not appear to scale 
well. Then, we can collapse the data in the scaling regime 
according to Eq. |(2Jl, merely using p* , v* , and /(oo) as 
free parameters to facilitate an accurate determination 
of y only. 

Fig. displays the same data as before, but only data 
judged by inspection to be suSiciently scahng {L > 1024 
at p = 0.35 down to i > 8 for p = 0.9, see Fig. H)) 
has been used for the fit. All data for p = 1 has been 
explicitly excluded due to the anomalous jump in the 
amplitude of a noted in Fig. 21 These choices are reflected 
in the collapse: All but p — 1-data combines exceedingly 
well for x » 1, predicting y = 0.256(1), only 1/2% above 
the exact value. For data scaling over almost 5 decades, 
one may have expected more accuracy for y, but slowly 
decaying scaling corrections [s^, and the smallness of y 
itself, limit the relative accuracy. The value for y remains 
quite robust under changes in the data points included 
in or excluded from the fit. On the other hand, the fit 
selects p* « 0.29, v* « 1.94, and /(oo) 2.1. Note that 
the fitted value for v* is larger than the actual value. It 
is this increase in v* that facilitates the collapse of the 
data in the asymptotic regime, x 3> 1. Consequently, 
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FIG. 7: Plot of the raw data for a{AE) as a function of 
systems size L for bond fractions p* < p < 1 using continu- 
ous Gaussian-distributed bonds. Some of this data is already 
shown in Fig. 1^1 
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FIG. 9: Collapse of the data from Fig. |7| according to the 
scaling Ansatz in Eq. (|5J. Here, the values for f{x), p* , v* , 
(j> and y are all determined from the fit of data in the scaling 
regime only. 
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FIG. 8: Collapse of the data from Fig. 

according to the 

scaling Ansatz in Eq. but using the exactly-known values 
for p* and v* (then y is not needed). Here, only /(oo) and (f> 
remain as free fitting parameters. 



the collapse of the data excluded from the fit at x < 1 
is somewhat poor (but still not too bad here). This is 
exactly the approach adopted in Refs. 113, ll^ to extract 
the best possible estimate for y. Hence, as noted there, 
the values for any parameter aside from y fitted in this 
way must be treated with caution. 

We can proceed in a similar manner for the bond- 
diluted hierarchical lattice with a continuous Gaussian 
bond distribution. First, we plot the raw data obtained 
with the algorithm described above in Fig.[7| In this case, 
finite-size corrections are more pronounced but appear to 
diminish more gradually towards p — 1. The scaling ar- 
guments from Sec, mil would suggest to fit the data to a 
form derived from Eq. 



/M 



(9) 



Here, it appears that the scaling collapse involves yet an- 
other parameter, the thermal-percolative cross-over ex- 



ponent (f> for the scaling variable x — L/S,{p) at p — > 
p* — Pc- In contrast to p* and v*, which can be deter- 
mined analytically for the hierarchical lattice, 4> is simi- 
larly nontrivial as y itself, and is even harder to estimate 
numerically. By definition, (j) — —vyp and yp has to be 
obtained exactly at p = p*, see Eq. (O, where the algo- 
rithm is the most delicate, limiting us to L < 2^^. As 
shown in Fig. O (bottom), from the data for a at p* we 
can extract yp — —0.9(1). With the percolation expo- 
nent, V = 1.22, we obtain (f) = 1.1(1). 

In Fig.|S| we collapse the data from Fig.[7|using Eq. Q 
by fitting / (oo) and 0, but holding p* and , and y fixed. 
Similar to Fig. |S1 the collapse proceeds well for data that 
is near x ^ 1 and has p near p* . Already for intermediate 
values of p, the data spreads widely, suggesting that ^(p) 
is too small there. We obtain a fitted value of </) « 0.9, 
not too far from the determination via yp above. 

Again, the focus on properties associated with p* gave 
us a good collapse near a; 1, but not about the desired 
asymptotic scaling regime for y at x ^ 1. Consequently, 
we again exclude all data in the fit that is by inspection 
of Fig. not yet in that asymptotic regime and proceed 
with an unrestricted fit involving all parameters. The re- 
sult of such a collapse is shown in Fig. |5| Now the data 
collapse is excellent in the asymptotic regime for data of 
sufficiently large L and large enough p, including p = 1, 
but poor near p = p* and for x < 1 for the excluded data. 
The fitted values here are p* « 0.32, u* « 1.20, cj) w 0.51, 
and y « 0.256, again within 1/2% error of the value di- 
rectly determined at p = 1. Surprisingly, even v* is well- 
approximated now while only (p is substantially off. In 
fact, ignoring the correction in the scaling behavior due 
to the exponent (f) (i. e. ^ = 0), the data collapse proceeds 
even more favorably, giving y = 0.2557 but p* — 0.31 and 
V* « 3. Having one less free parameter makes the data 
collapse a bit more robust, with v* picking up the error 
due to a less adequate scaling Ansatz. That a collapse of 
the data succeeds even when we ignore is a refiection 
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FIG. 10: Plot of the raw data for a{AE) as a function of 
systems size L for bond fractions pc — p* < p < 1 using 
continuous Gaussian-distributed bonds on a cubic lattice. At 
nearly each p, long transients are followed by an intermediate 
scaling regime, cut short at larger L, when numerical inaccu- 
racies result in a systematic drift. 



of the fact that in the limit of large argument for /, only 
two of the three exponents y, (f), and v* are independent. 
Eq. © reduces to cr ~ /(oo)L''(p-p*)'^+^''* , and setting 
= leads us right back to Eq. J^Jl. 



V. EDWARDS-ANDERSON MODEL WITH 
GAUSSIAN BONDS 

We have already pointed out the striking similarities 
between the data for defect energies on d = 3 hierarchical 
lattices here and on cubic lattices in Ref. for 
± J bonds. We complement this comparison here with a 
study of the EA on cubic lattices with Gaussian bonds. 
Such a study allows us to probe some of the assertions 
leading to Eq. . It also tests the universality of y with 
respect to the details of the bond distribution. The well- 
known result for d = I mentioned in the Introduction and 
recent studies for d > 1 |^ [s^l have show that there 
are significant differences in the results for the stiffness 
exponent with respect to bond distribution below the 
lower critical dimension, where j/ < 0, while Ref. [3^ 
has argued that universality should hold when y > 0. 
Our findings here, and for the hierarchical lattice above, 
support this point. 

Unlike for hierarchical lattices, the numerical effort 
require to achieve any reasonable system size L in the 
determination of defect energies on cubic lattices grows 
exponentially with the number of variables; the prob- 
lem is NP-hard |0. While the method of reducing low- 
connected variables is applicable just as well for the Gaus- 
sian distribution, our implementation of the extremal 
optimization heuristic does not perform well for these 
systems. Although it might be capable to yield reason- 
ably accurate predictions for ground states themselves, 
we notice a significant drift in the data for the defect 




L{p-p r 

FIG. 11: Collapse of the data from Fig. 1101 according to the 
scaling Ansatz in Eq. @. The values for the parameters had 
to be fixed by hand at p* — pc ~ 0.248, ~ 0.75, and 
S = 0.9. 




L(P-P )' 

FIG. 12: Collapse of the data from Fig. |7| according to the 
scaling Ansatz in Eq. i. e. (p = 0. Here, the values for 
f{oo), p* , u* , y are all determined from the fit of data in the 
scaling regime only, although all data from Fig. |7|are shown. 
Despite the obvious shortcomings of the data, the fitted value 
of y « 0.23 compares reasonably well with the best-known 
value of J/ = 0.24(1) [2^. l2^. drawn here for reference (line). 



energy, formed out of the difference between two closely 
related ground states, starting with remainder graphs of 
size n > 250. As Fig. 1101 shows, much of the obtained 
data, both for small and larger L, is not scaling and has 
to be discarded. We note that the long transients even 
before any asymptotic scaling is reached for small L is 
again reminiscent of the hierarchical lattice as in Fig. [3 
An attempt to collapse the data according to Eq. 
as before requires the additional effort in determining yp 
from Eq. A determination of the defect energy at 
Pc has the advantage that the exact reduction method 
almost always succeeds completely, obviating the appli- 
cation of any optimization and large lattice sizes can be 
reached j^^l- Our preliminary studies for systems up to 
L = 100 have yielded an exponent of yp = —1.28(2). 
The consensus of results for the correlation exponent in 
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d = 3 for percolation seems tohe v = 0.86(2) which 
results in a thermal-to-percolative crossover exponent, 
see Eq. (jHl, of (/) = 1.10(4). Note that while v < 1 
and \yp\ > 1 here, opposite to the hierarchical lattice 
in Sec. IIVI the values of (p are indistinguishable within 
errors. 

Similar to Fig. |S1 one may try to collapse the data 
according to the scaling Ansatz in Eq. lO by fixing p*, 
V* and y to their best-known values and fitting for f{oo) 
and 4>. Such a fit does not converge. We managed only 
to collapse the data for p ^ p* by hand with the best- 
known value for p* ^p^ra 0.248 [l^, but v* w 0.75, and 
(p = 0.9, somewhat below their best-known values. The 
result in Fig. ^2 while far less convincing for this limited 
data set, parallels that in Fig.|Hlto a large extend. 

Finally, wesel (j) — Q, which reverts Eq. Q into Eq. 
and eliminate all data that is obviously not scaling. Pur- 
suing a fit of the remaining data according to Eq. in 
the asymptotic regime, a; ^ 1, without fixing any pa- 
rameters. Fig. El is obtained. The result is noisy and 
apparently unsatisfactory: The value of « 3.1 is as 
large for Fig. 13 which separates this data in the scaling 
regime for each p into almost disconnected groups. Yet, 
the fit obtains a value of y « 0.23, in reasonable agree- 
ment with the best-known value oi y = 0.24(1) [27il2^. 

Clearly, as for the hierarchical lattice at the end of 
Sec. IIVI the focus on data with x 1^ i. e. p sufficiently 
larger than p* such that L :i> ^(p), makes the details of 
the transition aXp p* (such as the value of </>) irrelevant 
for the determination of y. This suggests Eq. as the 
appropriate scaling Ansatz to extract y, as proposed in 
Refs. [23,123, even for Gaussian bonds. 

VI. CONCLUSIONS 

We have explored a recently proposed method of 
collapsing data obtained on bond-diluted lattices to 
estimate low-temperature scaling properties for the 
Edwards- Anderson model. Using bond-diluted hierar- 
chical lattices from the Migdal-Kadanoff bond-moving 
scheme in d = 3, for which many properties at T = 
are known exactly or with high accuracy, the validity of 
the method is probed. The data obtained for the defect 
energy of the hierarchical lattice proves to exhibit the 
same qualitative features as the - much more limited - 
data for the EA, respectively for discrete and continuous 
bond disorder. 

The scaling Ansatz used to collapse the data, most 
generally Eq. (0), was proposed in the neighborhood of a 
T = phase transition at bond density p* between para- 
magnetic and spin glass behavior, which is closely associ- 
ated with the bond percolation transition at Pc |33l l3fi| | . 
Eq. (jSJ requires both, the divergence of the system size 



L and of the correlation length = {p — p*)~'^ ■ In 
contrast, the desired scaling behavior is associated with 
a r = fixed-point characteristic of the entire spin 
glass regime for L 3> 1, independent of ^(p) as long as 
L ^ £^{p). In particular, valuable data obtained for large 
L but intermediate values of p, for which ^{p) remains 
small, would have to be discarded in an Ansatz based 
on p — !■ p*. Thus, a naive Ansatz, Eq. (|3J), based on 
the limit L 3> £,{p) (i. e. x ^ 00) in Eq. exploits 
the obtained numerical data optimally (after data with 
x < 1 is cut). This Ansatz is "naive" in the sense that 
X = L(p — p* Y is not a true scaling variable, and the ob- 
tained values of the fitted parameters do not correspond 
to those defined by Eq. Those parameters provide 
useful degrees of freedom in the fit to remedy unknown 
corrections in ^(p), as our discussion here shows. In fact, 
Eq. Q has been used successfully in Refs. pT HM Iso f . 
where very robust scaling behavior was extracted for the 
desired defect scaling exponent y, but unrecognizable val- 
ues have been found for p* and v* . Unknown corrections 
in particular prevented a data collapse for data at p = 1 
there, a fact closely mirrored by the hierarchical lattice 
with a discrete bond distributions here, see Fig. ^ 

Clearly, the system sizes L attainable for the hier- 
archical lattice are unrealistic for the EA, and the un- 
known corrections to scaling disfavor any finite size scal- 
ing Ansatz in comparison to the accuracy obtained in a 
simple extrapolation of p = 1 data [35l lisf . Yet, the 
striking similarities in the behavior of the data for de- 
fect energies of the hierarchical lattice and EA should 
be noted. The Gaussian bond distribution results in ex- 
tended transient behavior until scaling is reached, making 
it impractical to extract the stiffness exponent y for the 
EA, even if an optimization heuristic would be applied 
that could handle remainder graphs beyond the limit of 
R:! 250 spins used here. Conversely, the cross-over at in- 
termediate p for system size corrections minimizes tran- 
sient behavior for ±J bonds, leading to the best data 
collapse for values of L realistic achievable for the EA. 
This intermediate window in p is invariant but narrows 
for d — > 00, allowing for reasonable determinations of yd 
for as high a dimension as d = 7, with scaling in the data 
collapse extending over two decades in d = 3 to merely 
half a decade in d = 7 . These new values for y^ allow 
for a direct comparison with mean field predictions |46j |. 
Our study here should add some confidence into those 
findings. 
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